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CN ■ ABSTRACT 



Aims. An optical search is carried out for supernova remnants (SNRs) in the Sc type nearby spiral galaxy M74, using ground- 
based observations at TUBITAK National Observatory (TUG, Antalya/Turkey) and Special Astrophysics Observatory (SAO, Russia). 
Observations are supplemented by the spectral analysis of archived X-ray data from XMM-Newton and Chandra. 
Methods. A survey of M74 covering ~ 9 arcmin- with [S II], Ho-, and their continuum filters. Interference filter images of M74 
are obtained with 1.5 m Russian Turkish Telescope (RTT150) at TUG and spectral data are taken with the 6 m Bolsoi Azimuthal 
Telescope (BTA) at SAO. The emission nebulae with continuum-subtracted line ratio values of [S II]/l/l67 16,6731 fHa > 0.4 are 

V^ ' ' identified as SNRs. A follow-up spectroscopy is obtained to confirm optical SNR identifications. 

^H ' Results. We have identified nine new SNR candidates in M74 with [S II]/Hq' > 0.4 as the basic criterion. We obtain [S III/Hq- ratio 

in the range from 0.40 to 0.91 and Ha intensities from 2.8 x 10"'^ erg cm"- s"' to 1.7 x lO"'"* erg cm"^ s"'. We also present spectral 
follow-up observations of these SNR candidates. However, we are able to spectrally confirm only three of them (SNR2, SNR3, and 

i-^ ' SNR5). The lack of confirmation for the rest might be due to the contamination by the nearby H II emission regions as well as due 

Mh' to the inaccurate positioning of the long slit on these objects. In addition, we search the XMM-Newton and Chandra Observatory 

archival data for the X-ray counterparts to the optically identified candidates. We find positional coincidence with only three SNR 
?** ! , candidates, SNRl, SNR2, and SNR8. The spectrum of SNR2 yields a shock temperature of 10.8 keV with an ionization timescale of 

■*^ ■ 1.6 X 10'" s cm"^^ indicating a relatively young remnant in an early Sedov phase which is not supported by our optical wavelength 

—^ ' analysis. Given the high luminosity of 10^' erg s"' and the characteristics of the X-ray spectrum, we favor an Ultra Luminous X-ray 

Source interpretation for this source associated with an SNR. We calculate an X-ray flux upper limit of 9.0 x 10"'^ erg cm"^ s"' for 
the rest of the SNRs including spectroscopically identified SNR3 and SNR5. 
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1. Introduction stars having a scale free nature in the galaxy disk. Fathi et al. 

,,_^ ,.,^- ,-„, . _ ■ , , ., ^o • (2007) have studied the internal kinematics of M74 over its en- 

M74 (NGC 628) is an Sc type spiral galaxy with a ~ 6 in- ^^.^ ^^^^ ^^. pabry-Perot interferometi-y with a good angular 

^ clination angle at an assumed distance of 7.3 Mpc (Shai'ina et resolution, and confirmed the presence of an inner rapidly rotat- 

O . !l-^^?^^' ^°"^ ^!: ^1;,.2004). It IS the brightest member of the sma ■ ^■^^_^.^^ component caused by the slow secular evolution of 

^ , M74 group in the Pisces constellation. M74 has spiral arms well ^j^^ large-scale spiral arms together with the oval structure. They 

- - formed with bright blue star clusters and dark cosmic dust lanes, ^j^^ ^^^^^^ ^^^^ 3qq j^ „ -^^^ -^ ^^^ ^^ ^-^^ calibrated 

In the early 2000, two supernovae (SN 2002ap and SN 2003gd) luminosities and diameters up to about 300 pc using continuum 

were identified in M74 and studied extensively. The supernova ^ubti-acted narrow band images in Ha. 

«J • event SN 2002ap is one of the four Type Ic SNe observed in X- 

d ■ rays (Van Dyk et al. 2003). SN 2003gd, which is found to be a ^74 is observed m the X-rays using the Chandra and 

■ ■ ' nearby Type II-P (plateau) event has a confirmed red giant pro- XMM-Newton data by Soria et al. (2004). 74 discrete X-ray 

eenitor (Hendrv et al 2005) sources are detected in the combined Chandra observations 

Leievre & Roy (2000) have recovered well over 100 small ^""J detection limi^t of * 6 x 10^^ erg s"'. They estimate 

and isolated H II regions and measured their fluxes in the ex- ^5% of the M74 discrete sources to have soft colors. They 

treme outer disk of M74 with R>R25 making use of deep Ha ^^f ^^I'^^l^t^ dermal SNRs with lurmnosities^ 2 x lO^J erg 

narrowband imaging. They show that massive star formation rate ll^T°r.^l'l^\ :17..T'''^' } ^J^^^^^'^^^^^^^' 

measured by the azimuthally averaged Ha surface brightness de- ^XOU J13651.1 + 154547) are detected by XMM-Newton m 

creases monotonically from the center out to R~20 kpc beyond ^^^ (Soria et al. 2004), while only the latter is detected by 

which it droDS raoidlv Chandra (Krauss et al. 2003). It is established that both of the 

Elmegreen et al. (2006) have examined the size and luminos- ^°^'"':^^ are variable. It is also argued that they can be Ultra 

ity distributions of the stai- forming regions in the galaxy with Luminous X-ray sources (ULX), however their true nature re- 

the Hubble Space Telescope. Their results suggest clumping of ^^™^ unclear 

Radio observations of M74 have been carried out during the 

Send offprint requests to: E. Sonbag Arecibo Galaxy Environment Survey (AGES). This survey is de- 
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signed to investigate M74 group environment in 21 cm with a 
better sensitivity and higher spatial and velocity resolution than 
previous observations (Auld et al. 2006). They obtain the spatial 
distribution of HI for M74 and other selected galaxies. 

In the UV regime, the surface brightness and color profiles 
for M74 are calculated by Cornett et al. (1994). They show that 
M74 disk has sustained significant star formation over the last 
five hundred million years. 

In this paper, we present our results on the search for new 
supernova remnant (SNR) candidates in M74 using the observa- 
tions by RTT150 and 6 m-BTA telescopes. We also make use 
of archived Chandra data to look for their X-ray counterparts. 
SNR studies are important in the theories of interstellar medium 
(ISM), and star formation since supernovae inject large amounts 
of matter, and energy into the ISM. However, in spite of the quite 
large number of Galactic SNRs their observations are impeded 
by several limitations such as the uncertainty in distances to indi- 
vidual objects and high extinctions along the line of sight in sev- 
eral regions of the Galactic plane. The limitations and uncertain- 
ties are inherently much less in extragalactic samples. Assuming 
that all SNRs are of the same distance to us for a given galaxy, 
we can easily compare their observed properties. Relative posi- 
tions of such SNRs are determined with more precision. Once 
we know the positions of SNRs, their distances with respect to 
H II regions and spiral arms can easily be calculated (Matonick 
& Fesen 1997; Blair & Long 1997, 2004). 

SNRs are identified in a number of nearby spiral galaxies 
using optical observations (e.g. D'Odorico, Dopita & Benveuti 
1980; Braun & Walterbos 1993; Magnier et al. 1995; Matonick 
& Fesen 1997; Matonick et al. 1997; Gordon et al. 1998,1999; 
Blair & Long 1997, 2004; Sonbas et al. 2009), X-ray observa- 
tions (Pence et al. 2001; Ghavamian et al. 2005), and radio ob- 
servations (Lacey et al. 1997; Lacey & Duric 2001; Hyman et al. 
2001). Multiwavelength surveys of SNRs are also carried out by 
Pannuti et al. (2000, 2002, 2007). 

In our work, we use a well known and accepted criterion (S 
llfHa > 0.4) proposed by Mathewson & Clarke (1973) to differ- 
entiate SNRs from typical H II regions. Rationale for this lies in 
the fact that in a typical H II region, sulfur is usually expected to 
be in the form of S ^^ due to strong photoionizing fluxes from 
central hot stars, making the ratio [S II]/Hq;, typically, in the 
range ~ 0.1 - 0.3. Shock waves produced by SN explosion prop- 
agate through the surrounding medium. The matter cools suffi- 
ciently behind these waves and variety of ionization states oc- 
cur including S^. This might be the reason for the increased [S 
II]/Hq' ratio observed in SNRs. It follows that almost all discrete 
emission nebulae satisfying the above criterion are expected to 
be shock-heated. 

The organization of the paper is as follows; In section 2, we 
discuss our imaging and spectroscopic observations and related 
data reduction. Search results in the optical band, the identifica- 
tion of SNRs, and search for their X-ray counterparts are pre- 
sented in section 3. Finally, conclusions and discussions of our 
results are provided in section 4. 

2. Observations and data reduction 

2.1. Imaging 

The imaging observations of M74 were performed with the 1 .5 
m Russian Turkish Telescope (RTT150) at TUBITAK National 
Observatory (TUG) in Turkey. M74 images were obtained using 
a Cassegrain imaging CCD Loral Lick 2k x 2k back-illuminated 
and anti-reflection coated (single layer, 700 A of hafnium ox- 



ide, Hf02), 2048 X 2048 pixel CCD with a plate scale of 0".26 
pixeL', giving a 9'.1 x 9'.1 field of view (FOV). We used nar- 
rowband interference filters centered on the lines of [S II], Ha, 
and continuum filters to remove starlight from Ha and [S II] im- 
ages. The interference filter characteristics used in these obser- 
vations are shown in Table 1 . Only one field image for each filter 
was obtained for M74 on September 19 and 20, 2004. We have 
exposure times of 12 x 300 sec for the [S II], Ha filters, and 4 x 
300 sec for the two continuum filters. This way, we were able to 
obtain deeper field images with higher signal-to-noise ratio for 
the faintest objects. 

The European Southern Observatory Munich Image Data 
Analysis System (ESO-MIDAS) software environment (version 
08FEBpll.l) was used to reduce the data (http://www.eso.org). 
We observed two or three standard stars each night from the list 
of Oke (1974) and Stone (1977) to determine the average flux 
conversion factors (FCF). We used the technique from Jacoby 
et al. (1987) to calculate this factor in the units of erg cm"^ s"'/ 
ADU s ' . FCF allows us to convert measured counts from an im- 
age to a flux in units of erg cm"^ s"'. Bias and flat-field frames 
for each image were also obtained. This way each exposure was 
bias subtracted, trimmed and flat-fielded. The cosmic hits were 
also removed from each [S II] and Ha image. Astrometry was 
applied to the resulting [S II] and Ha images to include world 
coordinate system information into the FITS header of the in- 
dividual images of M74. Red stars from USNO A2.0 catalog 
(Monet et al. 1998) were used to check positions of the stars. In 
Figure 1, we show the positions of SNR candidates overlaid on 
the Digitized Sky Survey (DSS) image of the M74 galaxy. 

2.2. Spectral data 

In order to resolve the [S II]/l/l6716,6731 lines and to separate 
[NII]/l/}6548,6583 lines from Ha providing a more accurate [S 
II]/Ha values than imaging ratio values, we carried out spec- 
tral follow-up observations. The spectral data of our SNR candi- 
dates were obtained by SCORPIO (The Spectral Camera with 
the Optical Reducer for Photometrical and Interferometrical 
Observations) which was mounted to 6m BTA with a 2048 x 
2048 pixel CCD (Afanasiev & Moiseev 2005). A 3500-7200 A 
spectral coverage at 10 A resolution and a 1" wide slit were used 
during the observations. 

The basic data reductions, flux and wavelength calibra- 
tions and interstellar extinction correction were performed 
using IDL codes and IRAF (Image Reduction and Analyses 
Facility) packages (http://iraf.noao.edu/docs/docmain.html). 
Spectrophotometric standard stars from Oke (1974) and Stone 
(1977) catalogs were also observed each night. We derived 
fluxes from the spectrophotometric standards and used these 
to calibrate the fluxes for spectral lines in our SNR spectra. 
Biases, internal lamp flats, and Ne-Ar calibration lamp frames 
were obtained for each observation set. The splot IRAF routine 
was used to measure the emission line fluxes of [S II] and Ha. 
While obtaining spectral data of SNR candidates, a bright star 
was placed on the slit along with the target object. Galactic 
extinction correction was applied to the SNR spectra as part of 
the spectral analysis (Cardelli et al. 1989). 

3. Results and discussion 

3. 1 . Results from Optical Imaging Analysis 

We used the technique described by Blair & Long (1997) in 
order to identify the SNR candidates. The identification tech- 
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nique was constructed through the comparison of continuum - 
subtracted [S II] /1/167 16,6731 and Ha images. Eventually, emis- 
sion nebulae that have region-integrated values of [S II]/HQf> 
0.4 were identified as the SNR candidates. Preliminary search 
to find these candidates were carried out through a comparison 
of the continuum subtracted [S II] and the Ho- subfield images 
by a blinking technique. In the process, a full-field image of the 
galaxy was divided into regions of about 2' by 2' squares for the 
visual inspection and assessment of the fields for the candidates. 

When a point like source satisfied the criterion [S II]/!!!]- > 
0.4 we then marked it as an SNR candidate. After the assessment 
of all 2' by 2' subfields of M74, continuum subtracted images 
were used to determine the ratio of [S II]/!!^ counts. During this 
process, background counts were subtracted from the [S II]/Ha 
ratio images using annuli farther away, but centered on the se- 
lected candidates. The background extraction areas were also 
normalized to the source extraction areas. At the assumed dis- 
tance of 7.3 Mpc for M74, our aperture size corresponds to ~ 
50 pc. During our interference filter imaging observations see- 
ing was 1.4", thus we can not extract any radius information for 
our SNR candidates since it was limited with the seeing for the 
night. 

We found nine emission nebulae consistent with [S II]/Ha 
> 0.4 criterion for SNR identification with the technique de- 
scribed above. The list of new SNR candidates with corrected 
flux ratios are given in Table 2. Ha fluxes for SNR2, SNR3, and 
SNR5 are taken from spectral measurements while remaining 
ones are from image data. A limiting flux sensitivity of ~ 10"'^ 
erg cm~^s"' was calculated by choosing a structure with a limit- 
ing magnitude in the field of the galaxy. M74 covers a total area 
of about lO'.S X 9'.5 in the sky. Since FOV of our CCD is 9'.1 
X 9'.1, each filter image covered almost the full face of M74. 
Therefore only a single observation for each filter was required. 
We zoomed on the full continuum-subtracted Ha image of M74, 
and extracted four subfields, each of size 3'. 5 x 2'. 5, and marked 
our SNR candidates on these images as displayed in Figure 2. 

3.2. The Follow-up Optical Spectroscopy 

We also conducted spectral observations to confirm the optical 
SNR candidates detected by optical imaging analysis. Among 
the nine candidates detected in M74, we were able to confirm 
only three of them (SNR2, SNR3, and SNR5) with the derived 
specific line ratios [S II]/!!^ of 0.46, 0.91, and 0.58 respectively. 
Their optical spectra are shown in Figures 3-5. Observed line in- 
tensities relative to H/3, the E^b-v) values and Ha intensities for 
the spectra of these three SNRs are given in Table 3. The rest 
of the six SNR candidates (SNRl, SNR4, SNR6, SNR7, SNRS, 
SNR9) could not be confirmed using the spectroscopic line ratio 
values. Several factors could be affecting the spectral [S II]/Ha 
line ratio values. For instance, the long slit used for light acquisi- 
tion may not always be accurately pointed on the objects. In such 
cases contamination from nearby or overlapping H II emission 
regions could be affecting the line ratios. This seems likely for 
SNRl, SNR6, SNR7, SNRS, and SNR9 that were quite close to 
extensive H II emission regions as shown in Figure 6. The spatial 
coordinates and sizes of these H II regions were taken from the 
catalog of Fathi et al. (2007) (The catalog can also be reached 
at http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i-A/466/905/). We 
thought that there could be significant H II emission contamina- 
tion into these sources, as was foreseen by Matonick & Fesen 
(1997), and Blair & Long (2004). They discuss some of their 
SNR candidates in relation to closeby H II regions (e.g., their 
SNR5 in NGC6946, SNR39 in MlOl, and SNR19 in MSI) with 



the angular separations of the SNRs being up to 15" - 40" from 
H II regions. The angular separations of our candidates SNRl, 
SNR6, SNR7, SNRS, and SNR9 from the nearby H II regions 
calculated using our images were found in the range 5" to 10" in 
M74. This corresponds to approximately the same distance scale 
when the relative distances of target galaxies are considered. On 
the other hand, we estimate that SNRl, SNR4, SNRS, SNR9 
have significant total [O I] emissions reaching about from 20% 
to 60% of the H/3 intensity. This is considered to be a secondary 
shock-heated gas indicator. In turn, presence of a detectable [O 
I]/16300 emission is considered to be a viable argument for SNR 
identifications, in general (Blair & Long 2004). 

We note that, the detected SNR distribution in M74 appears 
mostly in the southern field of the galaxy. Similar asymmetric 
distributions are observed in some other spiral galaxies such as 
NGC 2403 and NGC 2903 (Matonick et al. 1997; Sonbas et 
al. 2009). Presently, the leading explanation for such asymme- 
tries is that in the half of a given galaxy where H II regions are 
brighter, less SNRs are detected. This seems to be the case also 
for M74. We believe that this is not an artifact of our imaging 
or SNR identification techniques. We propose that, such asym- 
metries need to be discussed further with more galaxy examples 
and observations. 

At the most basic level, average electron density in a nebula 
is usually measured through effects of collisional de-excitation. 
This is usually done by the comparison of intensities of two lines 
from the same ion, emitted by different levels with the same 
excitation energy. One of the best examples of this procedure 
makes use of the electron density ratios [S II]/16716/[S II]/16731 
and [OII]/13729/[OII]^3726. In our work we used [S II]^6716/[S 
II]/16731 line ratio from SNR spectra to calculate the electron 
density, Ng . The Space Telescope Science Data Analysis System 
(STSDAS) task nebular.temden, based on the five-level atom ap- 
proximation, calculates the required electron density given the 
electron temperature. Our calculated A^^ value is 190 + 100 cm"-' 
for SNR3 corresponding to the line ratio of 1.24 for an assumed 
electron temperature of T = 10'* K. The reason for the error term 
quoted here arises from the fact that the theoretical model loses 
its sensitivity for electron densities < 100 cm"^ which is in line 
with our spectral measurement sensitivities. With the same rea- 
soning, although the Hne ratio values of [S II]/16716/[S II]/i6731 
> 1.46 is considered to correspond to a low density limit with 
Ne < \Q cmT^ (Osterbrock 19S9), we can only quote an elec- 
tron density for both SNR2 and SNR5 of simply < 100 cm ^ 
Using the shock abundance models by Dopita et al. (19S4), the 
elemental abundances can be determined from the SNR spec- 
tra. Their models assume a shock velocity around or less than 
100 km s ' which is considered a low shock velocity. Same au- 
thors also make use of the fact that the behavior of [OIII] line is 
the most sensitive indicator of shock velocity variations. These 
variations fall rapidly for values less than S5 km s"' . In our case, 
[OIII] ^5007/H/? emission Hne ratios for SNR2, SNR3 and SNR5 
were found as 0.27, 1.22 and 1.6, respectively, corresponding 
to nebular propagation velocities of < S5 km s"' (see Figure 5 
of Dopita et al. 19S4). There are many SNRs which have low 
shock velocity values in a number of nearby galaxies (Matonick 
& Fesen 1997; Blair & Long 2004). 

Among the working methods to separate SNRs from H II re- 
gions are the plots of [S II]/Hq; versus H/3 flux and [S II]/Hq; 
versus [NII]/Hq' line ratios. It was noted that in a number of 
nearby galaxies many SNR candidates show enhanced [Nil] and 
very weak H/3 emissions when compared with the H II regions. 
We adopted the method by Blair & Long (2004) based on their 
Figure IS where the same two plots for SNRs in MS3 are given. 
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When we inserted our ratio values for detected SNRs into these 
plots (Figure 7), we find similar tendencies. Our new SNR can- 
didates are also located in the lower end as would be expected. 
The observed low li/3 flux values are indicative of faint ionized 
hydrogen regions. 

3.3. Comparison with X-ray Observations 

In a cross-check with the X-ray point source list of 
M74, three X-ray sources (CXOUJ013631.7+154848, 
CXOUJ013644.0+154908, CXOUJ013646.0+154422, Soria 
et al. 2004) are found coincident with our optically detected 
SNRs (SNRl, SNR2, SNRS) taking ~ 8" positional error circle 
around the objects. The error in association of SNRl and SNRS 
is ~ 5" and SNR2 is ~ 3" . In Table 4 we list the optical and 
X-ray associations for these SNRs. Among these, only SNR2 is 
a spectrally confirmed candidate. We used an archival Chandra 
pointed observation of M74 obtained in 2001 June 19 (with an 
exposure of 46.3 ksec) to derive a spectrum of the SNR. We 
performed standard point source spectral analysis on the data 
and extracted source photons from the location of SNR2 with 
a circular extraction region of 12". 5 and background photons 
were extracted from a region normalized to the source photon 
extraction region. Entire analysis was performed using CIAO 
3.4 (http://asc.harvard.edu/ciao34/indexJitml) and XSPEC 
12.5.0ac software (Arnaud 1996). 

We measured a count rate of (1.5±0.5) x 10""' cts s"' for 
SNR2. We derived a source spectrum and fitted this spectrum 
using the XSPEC model PSHOCK, a non-equilibrium ioniza- 
tion plasma emission model for shock-heated plasma appropri- 
ate for SNRs (Borkowski et al. 1996). The final Chandra spec- 
trum of SNR2 is displayed in Figure S. The resulting spectral 
parameters are Nh=1.1 J""* x 10^^ cm^^, kTp,as,„a=l0.^1f;f keV, 
the ionization timescale t =not=1.6>''^ x 10'" s cm"^, and the 
model normalization K=l.l>^^ x 10"^ cm"^; the lower limit 
of the normalization is unconstrained. The;^^^ of the fit is 1.3. 
The unabsorbed source flux is about 1.7 x lO"'-' erg cm"^ s"' 
which translates to an X-ray luminosity ~ 1.1 x lO-'^ erg s"' at 
the source distance of 7.3 Mpc. The derived plasma temperature 
yields a shock speed of 3 100 km s"' (kTs = (3/ 16)//m//(vj)^) us- 
ing Rankine-Huguenot jump conditions. The model normaliza- 
tion can be used to calculate the electron density of the shocked 
material as n(.=0.45 cm"-' (filling factor f=l) given a volume lim- 
ited by the seeing that corresponds to 50 pc in size at the source 
distance. This value is consistent with the density derived from 
our optical line ratio analysis n^, < 10 cm"-'. Assuming that the 
ambient density no is not going to be more than the electron den- 
sity, the ionization timescale yields a minimum life time of about 
1120 years for SNR2. This indicates that SNR2 maybe in the 
Sedov phase of its evolution. 

It is important to note that the calculated luminosity for 
SNR2 is large for an SNR in the Sedov phase of evolution and 
the calculated expansion speed can not be supported with the 
blue shifts derived from the optical spectroscopy. Given this re- 
sult, we also suspect that SNR2 may be a ULX (e.g. XMMU 
J013636.5+155036; Soria & Kong 2002). We modeled the spec- 
trum of SNR2 with a power-law yielding an N//=2.6J"-^ x 10^' 
cm"^ a r=1.8+^-^ keV and a normahzation Kpi^3.lf'' x lO"*^ 
with axl-^.^. Adding a blackbody component along with the 
power-law gives a kT ~ 0.6-2 keV and steepens the power-law 
index to l.S-5.4. A single blackbody fit gives a large ;t'y above 
1.5. 



Given the high shock plasma temperature yielding expan- 
sion speed that can not be supported with the optical data and 
that the X-ray luminosity is of the order of 10^^ erg s"', we do 
not favor an old SNR interpretation for SNR2 as expected from 
the optical data analysis. We favor a ULX interpretation associ- 
ated with a SNR. We expect that the very slow (less than S5 km 
s"') expansion speed from the interpretation of the optical data 
will not yield detectable X-ray emission at the distance of M74. 
Then, any emission resulting from the point source remnant of 
the SN explosion, will dominate the X-ray emission, though the 
SNR may dominate the optical emission. The ULX is expected 
to rejuvenate and ionize the nebulosity around it. Associations 
between SNRs and ULXs have been found or suggested before 
(Roberts et al. 2003, Abolmasov et al. 200S, Di Stefano & Kong 
2004; Di Stefano, et al. 2010). 

We also checked the other two sources, SNRl and SNRS. 
These sources were detected with (1.0+0.5) x 10""* cts s"' and 
(2.7+0.5) X 10"^ cts s"', respectively. The fits with the PSHOCK 
model the two sources result in a neutral hydrogen absorption 
of about 0.9 X lO^^ cm"2, and a kTp/„s„,„ of about 34-73 keV 
with an ionization timescale r in a range 2-5 x 10'" s cm"-'. The 
spectral parameters yield an unabsorbed X-ray flux of 5 x 10"'^ 
erg cm"^ s"' and 1 x 10"'^^ erg cm"^ s"' for SNRl and SNRS, 
respectively, translating to an X-ray luminosity of ~ 5 x 10"'^ 
erg s"' and ~ 7 x 10-'** erg s"'. Given these spectral parameters 
and assuming these are SNRs we can speculate that the expan- 
sion speed is 2-3 times larger than SNR2 for these two SNRs, 
the electron density is about 1.4 to 2.3 times less than SNR2, 
and the minimum life time is about 1570-2600 years for SNRl 
and SNRS compared with SNR2. In general, one expects that 
the SNR shock decelerates in time, thus these longer life times 
are in contradiction with the faster shock speeds and high X-ray 
temperatures detected for these two SNRs in comparison with 
SNR2. These shock temperatures and expansion speeds are not 
supported by our optical data. We note that these two sources 
are not confirmed SNRs with our ratio measurements using op- 
tical spectroscopy. Thus, we suspect them as SNRs, we favor an 
X-ray binary or an AGN interpretation for the two sources. For 
the same three sources XMM-Newton data yielded confirming 
spectra. Therefore, no details are repeated here. 

Finally, we calculate an X-ray flux limit of < 9.0 x 10"'^ erg 
cm"^ s"' with a luminosity limit of < 5.0 x 10^^ erg s"' for the 
rest of the SNRs in our list including the two other spectrally 
confirmed sources SNR3 and SNR5 since we find no positional 
coincidence with any of the source listings of M74. 



4. Conclusions 

We conducted a survey of SNRs using optical imaging and spec- 
troscopic measurements in the nearby spiral galaxy M74. In this 
survey, we used blinking between continuum-subtracted Ho- and 
continuum-subtracted [S II] images to deduce SNR candidates 
for a given ratio between the lines. The SNRs were confirmed 
by additional spectroscopic observations and finally, comparing 
the optically detected SNRs with archived Chandra observations. 
Nine SNRs were identified in M74 by our method. Three of these 
have shown positional coincidences with X-ray sources. SNR2 is 
our best candidate to be an SNR using optical wavelength anal- 
ysis and the X-ray source is consistent with a ULX interpreta- 
tion, making it one of the few associations between a SNR and 
a ULX. We calculated an X-ray flux upper limit of 9.0 x 10"'^ 
erg cm"^ s"' and a luminosity upper limit of 5.0 x 10^^ erg s"' 
for our new SNR candidates excluding SNRl, SNR2 and SNRS. 
Our spectroscopically detected SNRs have shock velocities < S5 
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kms ' indicating that they fall into a low shock velocity range 
implying very old remnants from which little X-ray emission is 
expected. 

According to the SN rate calculation of Matonick & Fesen 
(1997) in such galaxies, about half of all SNe are of Type Ib/c 
or Type II, produced by massive stars. In turn, only half of these 
SNe are located in regions with enough ambient density to pro- 
duce a detectable SNR. This means that, only about a quarter of 
all SN events may leave easily detectable optical remnants. In 
this respect four times more SN events should have exploded in 
M74 and we detected nine SNRs. With this reasoning, we would 
expect to see ~ 40 SNRs. If the optically observable lifetime of 
an SNR is assumed to be ~ 20,000 years (Braun, Goss, & Lyne 
1989), a crude SN occurrence rate for M74 can be estimated giv- 
ing a value of about 1 per ~ 500 years. In the light of observation 
of two recent SNe from M74, our result can be reconciled if we 
accept the fact that more SNe occur at places that are difficult to 
reach optically, (i.e. in gas-rich regions of the galaxy). 
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Table 1. Characteristics of the interference filters used in our 
observations 



Name 


A 


FWHM 




Wavelength (A) 


A 


[SII] 


6728 


54 


continuum 


6964 


350 


Ha 


6563 


80 


continuum 


6446 


123 



Table 2. New optical SNR candidates detected in M74 in our 
observations. 



SNR Name 


RA 


DEC 


[S Ill/Ha 


I(Ha) 




(J2000.0) 


(J2000.0) 




(erg cm ^s ') 


SNRl 


01:36:43.8 


+ 15:49:12.7 


0.47 


1.5E-14 


SNR2 


01:36:31.5 


+ 15:48:50.3 


0.46 


6.0E-15 


SNR3 


01:36:46.9 


+ 15:45:15.2 


0.91 


2.9E-15 


SNR4 


01:36:45.3 


+ 15:43:06.6 


0.49 


1.4E-14 


SNR5 


01:36:48.9 


+ 15:43:34.4 


0.58 


2.8E-15 


SNR6 


01:36:45.1 


+ 15:45:25.3 


0.47 


1.5E-14 


SNR7 


01:36:45.9 


+ 15:44:46.1 


0.45 


1.6E-14 


SNR8 


01:36:45.9 


+ 15:44:27.1 


0.46 


1.6E-14 


SNR9 


01:36:42.4 


+ 15:44:16.3 


0.40 


1.7E-14 



'i»*»l(rt 




[oin] 



HP 



iKh^ 



[MI] 



[HI] 
\ 



^^X^u 



[SII[67I6 



/ 



[Sn]6731 



w*! 



4000 5000 6000 7000 

l(WaYeIeii^tli).A 

Fig. 3. Optical spectrum of SNR3 in M74 obtained with BTA. 
Identified lines are indicated. 
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Fig. 4. Optical spectrum of SNR5 in M74 obtained with BTA. 
Identified lines are indicated. (The strong spike just blueward of 
Hy is most probably due to an unsubtracted cosmic hit on the 
CCD chip in use.) 
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Fig. 5. Optical spectrum of SNR2 in M74 obtained with BTA. 
Identified lines are indicated. 
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Fig. 8. X-ray spectrum of SNR2 obtained using archival Chandra 
data in 0.3-10.0 keV fitted with a neutral hydrogen absorp- 
tion model, WABS and a PSHOCK non-equilibrium ionization 
plasma emission model (WABS'-PHSOCK). The crosses indicate 
the data, the solid line is the fitted model(upper panel) and the 
residuals between the data and the model in standard deviations 
(lower panel). 



E. Sonbag et al.: A Search For Supernova Remnants in The Nearby Spiral Galaxy M74 (NGC 628) 



5NR1 
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Fig. 1. Nine new SNR candidates found in this work are labeled with circles on the image of M74 extracted from Digital Sky Survey 
(DSS). 





Fig. 2. Positions of SNRl (a), SNR2 (b), SNR3, SNR6, SNR7, SNR8, SNR9 (c), and SNR4, SNR5 (d) overlaid on a ~ 3 '.5 x 2 '.5 
subfields of continuum-subtracted Ha images of M74. 
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Table 3. Relative line intensities and observational parameters for SNRs spectroscopically observed in M74 



Line 


SNR2 


SNR3 


SNR5 


H/3(/14861) 


100 


100 


100 


0111(^4959) 


3 


60 


50 


0111(^5007) 


27 


122 


160 


NII(^5200) 


- 


3 


- 


He(i5876) 


13 


33 


9 


OI(-16300) 


17 


21 


32 


OI(-16364) 


14 


- 


7 


NII(^6548) 


30 


154 


56 


Hq'(^6563) 


288 


589 


335 


NII(^6583) 


48 


281 


85 


SII(^6716) 


92 


297 


122 


SII(^6731) 


41 


239 


73 


E(B-V) 


0.07 


0.06 


0.06 


I(Ha) 


6.0E-15ergcm-2s-i 


2.9E-15ergcm-2s-i 


2.8E-15ergcm-2s-i 


[SII]/Ha 


0.46 


0.91 


0.58 



Table 4. Spatial associations between optically identified SNRs 
and Chandra-detected X-Ray sources (from the combined 2001 
June - October observation). 

Optically detected SNR RA DEC X - Ray Source 

SNRi 01:36:43.8 +15:49:12.7 CXOUJO 1 3644.0+ 1 54908 

SNR2 01:36:31.5 +15:48:50.3 CXOUJ013631.7+154848 

SNR8 01:36:45.9 +15:44:27.1 CXOUJ013646.0+154422 
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Fig. 6. H II regions and detected SNRs of M74. Perpendicular short lines indicate SNRl, SNR2 and SNR3 that are associated 
with ChandraX- ray point sources CXOUJ013631.7+154848,CXOUJ013644.0+154908,CXOUJ013646.0+154422, respectively. 
White circles indicate H II regions, and black circles indicate our six other SNRs. 
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Fig. 7. [S II]/Hq' vs [NII]/Hq' ratios and [S II]/Ha' vs HjS fluxes. Stars indicate our spectroscopically identified SNRs while diamonds 
those unidentified ones. These are compared with SNRs of M83 Blair & Long (2004) where open circles refer to spectroscopically 
identified SNRs, semifilled circles refer to spectroscopically unidentified SNRs, and filled circles refer to H II regions. Similar 
tendencies between SNRs of M83 and M74 are noticeable. 



